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We derive a hydrodynamic model for a liquid of arbitrarily curved flux-lines and vortex loops 
using the mapping of the vortex liquid onto a liquid of relativistic charged quantum bosons in 2+1 
dimensions recently suggested by Tesanovic and by Sudb0 and collaborators. The loops in the flux- 
line system correspond to particle-antiparticle fluctuations in the bosons. We explicitly incorporate 
the externally applied magnetic field which in the boson model corresponds to a chemical potential 
associated with the conserved charge density of the bosons. We propose this model as a convenient 
and physically appealing starting point for studying the properties of the vortex liquid. 



I. INTRODUCTION 



The physics of vortex matter has been a very active 
field of research since the discovery of high-temperature 
superconductivity Jl|,||. The melting of the Abrikosov 
vortex lattice at a field well below the mean field upper 
critical line, H C 2(T), is now well established both theo- 
retically and experimentally. The nature of the resulting 
vortex liquid phase remains, however, unclear. Various 
approaches have been used to study the properties of 
the liquid, ranging from the mapping of the statistical 
mechanics of vortex lines onto that of (nonrelativistic) 
two-dimensional quantum bosons |][(| , to continuum hy- 
drodynamic models . Recently, it has been proposed 
that the vortex liquid phase should be viewed not as a 
collection of well-defined directed vortex lines induced by 
the externally applied field, but rather as a system where 
closed vortex loops proliferate, with components in all di- 
rections p!p| , p3| . This picture has been substantiated by 
numerical work by Sudb0 and collaborators Jl2|-p8[ . Such 
a vortex liquid has been described by Kiometzis et al. JTg| ] 
and by Tesanovic pOj , pd[ | by mapping it onto a system of 
relativistic two-dimensional quantum bosons, where the 
proliferation of vortex loops and overhangs in the flux- 
line system corresponds to particle-antiparticle creation 
and annihilation events in the bosons. This mapping has 
been developed in some detail for the case of zero external 
field. 

In this paper we consider a vortex-line liquid in an ar- 
bitrary homogeneous external field, H. The only restric- 
tion is that H be well below the mean field upper criti- 
cal line, H C 2, so that the London approximation applies. 
The liquid in general consists of both field-induced vor- 
tices (on the average aligned with the external field, but 
where large fluctuations leading to overhangs are not ex- 
cluded) and closed vortex loops generated spontaneously 
by thermal fluctuations. On the basis of general consid- 
erations, the statistical mechanics of such a vortex liquid 



maps onto that of relativistic charged quantum bosons in 
2+1 dimensions coupled to a screened electromagnetic 
field. Our model generalizes that of Tesanovic Jll]] in 
that it explicitly incorporates the external field H which 
in the boson model corresponds to a chemical potential, 
fj, r , associated with the conserved charge density of the 
bosons. By manipulating the relativistic boson action, 
we then derive a hydrodynamic model for a liquid of ar- 
bitrarily curved flux- lines and closed vortex loops. We 
propose this model as a convenient and physically ap- 
pealing starting point for studying the properties of the 
vortex liquid. 

We consider for simplicity a clean, isotropic type-II 
superconductor in the mixed state. The anisotropic case 
is discussed in the Appendix. In the London limit, the 
familiar Ginzburg-Landau free energy functional can be 
rewritten in terms of discrete vortex degrees of freedom 
by parametrizing the z-th line by its position i*i(sj), with 
Si the arclength along its curve, as j2(^-^2j 

g[{r,( Sl )},H] = Uself + U lnt 
a N r 

i=l J 

where 

JV 

E/^/weiJ] / dsi (1.2) 

2—1 

is an approximate form for the self-energy of a single vor- 
tex line, and 




with V{r) — (eo/r)e~ r / x , describes the screened inter- 
action between vortex segments on different lines. Here 
€q = 0§/(47rA) 2 and ei = eglnK. Also, A is the effec- 
tive London penetration depth and k = A/£, with £ the 



*Present address: Physics Department, Harvard University, Cambridge, MA 02138 



1 



coherence length. The prime on the integral sign indi- 
cates that the interaction must be cutoff at intervortex 



separations of the order of £. The last term in Eq. (1.1) 
represents the interaction with the homogeneous external 
field, H. A short distance cutoff ks £ has been assumed 
to regularize the integrals at short scales. The canonical 
partition function of the system is 



Z(H,T) = J V{ 



Ti(Si)} 



-Q/k B T 



(1.4) 



The prime over the integral sign indicates that the func- 
tional integration must be performed over all flux-line 
configurations {vi(si)} subject to the constraint V • B = 
(i.e. no flux-line can start or stop within the sam- 
ple). The expression for the vortex energy given in Eqs. 
(1.1-1.3) is the starting point of much theoretical work 
on vortex arrays. As required, the free energy is a scalar 
and it is therefore rotationally invariant. 

A familiar approximation to the general free energy 
given in Eq. (1.1) can be obtained at low fields and 
temperatures, where field-induced vortex lines interact 
weakly and fluctuations away from the external field 
direction are small. In this case it is convenient to 
choose the z axis in the direction of the applied field, 
H = zHq and to parametrize the vortex positions as 
r, = [rj^(z), z] . Asuming small transverse fluctuations 
of the vortices and retaining only the interaction among 
vortex segments at the same "height" z (to be referred 
to as the local approximation), one obtains 



N 



g[{ ri ( Sl )},^o] -Y, dz 



dr ±i (z) 



dz 



i¥=3 

with n = #000/(471-) - ei HJE). Eq. (fj) is clearly only 
appropriate to describe field-induced directed lines and it 
neglects the possibility of thermally-induced closed vor- 
tex loops. It is well known that the partition function of 
the array of directed vortex lines described by Eq. (1.5) 
maps onto the quantum mechanical partition function in 
the path integral formulation for a gas of two-dimensional 
bosons H|]j2^-|2j| . In this mapping, the vortex lines cor- 
respond to the boson world lines, e± is the boson mass, 
and z represents the imaginary time. This boson map- 
ping has been exploited by Nelson and collaborators to 
study the properties of vortex liquids and is particularly 
useful for the description of vortex arrays pinned by cor- 
related disord er p ] . The approximate vortex free energy 
given in Eq. (1.5) is invariant under a uniform tilt of the 
external field away from the z direction. In the boson 
analogy, this translates into invariance with respect to 
Galileian boosts |^6| . 

One of the limitations of this boson model is that it ne- 
glects the nonlocality of the intervortex interaction in the 



external field direction. Feigel'man and collaborators ar- 
gued that this nonlocality can be incorporated in the bo- 
son formalism by mapping the vortex partition function 
onto the path integral of two-dimensional charged bosons 
p7[ coupled to a screened gauge field that mediates the 
non-instantaneous interaction among the bosons p8| , p9[ . 
The boson action proposed by these authors suffers, how- 
ever, from one crucial problem — it does not incorporate 
the symmetry of the original vortex free energy given in 
Eq. (1.1). More precisely, the boson analogue of a ro- 



tation of the three-dimensional space where the vortices 
are embedded is a Lorentz transformation of the 2 + 1- 
dimensional space-time where the bosons are embedded. 
Since the original vortex energy is rotationally invariant, 
its dual boson action should be Lorentz invariant. But 
the model implemented by Feigel'man and collaborators 
lacks a definite symmetry as the field and the interaction 
part of the Lagrangian density are Lorentz invariant, but 
the approximate form used for the free-particle part is 
only Galileian invariant. We note that the non-Gaussian 
hydrodynamic model proposed recently by us || suffers 
from the same problem, as it was actually derived from 
the nonrelativistic charged boson model. This problem 
is in retrospect immediately apparent in the hydrody- 
namic formulation, where the free energy can be written 
in terms of the coarse-grained fields [h H (r) , t H (r)l as the 
sum of three terms, as in Eq. fll.ip , 



Q H [n H ,i H ,H ] = U« lf + U& - NLH Q ^- , (1.6) 
with 



h 1 f |t H (r)| 2 



2j r 1 n H (r) 



(1.7) 



and 



U H = — 



K t {v-r')i H {v)-t H {v') 
■f K c {r-v')n H {Y)n H {v') . (1. 



Here K c (v) is the real space compressional modulus and 
K t (r) is the compressional part of the real space tilt mod- 
ulus. For an isotropic superconductor, K c (r) — K t (r). 
Since [n H (r),t H (r)] is a vector under rotations (pro- 
portial to the local induction B in the superconduc- 
tor, with n H (r) ~ B z (r) and t H (r) ~ Bj_(r)), we see 
that, for L — ► oo, U^ t is a scalar and therefore rota- 
tionally invariant, while is not a scalar and there- 
fore is not rotationally invariant. Note that, in con- 
trast, both the uncharged boson model a la Nelson |25j ] 
and conventional Gaussian hydrodynamics p(i|,f7t (where 
the field n H (r) in the denominator of Eq. (1.7) is re- 
placed by its equilibrium value) are consistent approxi- 
mations to the free energy of directed vortex lines and 
respect the restricted symmetry under global rotations 



2 



about the z axis. These observations indicate that, as 
recently pointed out in Refs. ]l9| , [To| - |T5| , ^8t , the proper 
two-dimensional boson model dual to the original three- 
dimensional vortex problem is a gas of relativistic charged 
bosons in 2 + 1-dimcnsional space-time, interacting via a 
screened electromagnetic field. 

In this paper, starting from the action of relativistic 
charged quantum bosons, we derive a nonlocal, nonlin- 
ear hydrodynamic model of liquids of arbitrarly curved 
vortex lines that properly incorporates the full rotational 
symmetry of the vortex energy. As discussed in Refs. 
ju^lO - 18 1 , this leads to some interesting new concepts, 
namely, the appearence of thermally excited antivortex 
lines and vortex loops. It is well known that a loop gas 
is dual to a gas of relativistic bosons in 2 + 1 space-time 
dimensions. This duality has been exploited recently by 
a few authors to study the properties of vortex liquids 
in type-II superconductors. Kleinert and collaborators 
p9| , PH have developed a field theory of a gas of fluc- 
tuating closed vortex loops and have applied renormal- 
ization group ideas to study its phases. Tesanovic has 
proposed a theory |ll]] which separates "field-induced" 
from "thermally-generated" degrees of freedom. The for- 
mer correspond to directed flux-lines, while the latter 
correspond to closed vortex loops. A prediction of this 
work is the existence of a phase transition, dubbed the 
^-transition, within the flux liquid phase from a low 
temperature liquid of lines with finite tension to a liq- 
uid phase that lacks superconducting coherence in all di- 
rections and where vortex loops proliferate yielding the 
vanishing of the flux-line tension. Numerical evidence for 
this ^-transition has recently come from the numerical 
simulations of Sudb0 and collaborators |Q . 

In this paper, we propose a long wavelength descrip- 
tion of a liquid of interacting field-induced vortex lines 
and vortex loops which treats both on equal footing. The 
external magnetic field is incorporated explicitly and con- 
trols the density of field-induced lines. In the dual de- 
scription it enters as the chemical potential of the bosons, 
in analogy to the nonrelativistic boson model. The im- 
portant difference here is that the chemical potential cou- 
ples to the local charge density of bosons, which is the 
relevant conserved variable as contrasted to the parti- 
cle number density. The hydrodynamic model proposed 
here is physically appealing model and naturally gener- 
alizes the familiar Gaussian hydrodynamics. Its ramifi- 
cations and particularly its standing with respect to the 
aforementioned ^-transition are interesting directions for 
future work. 



II. DERIVATION OF THE HYDRODYNAMICS 
OF FLUX-LINES AND FLUX-LOOPS FROM A 
RELATIVISTIC BOSON ANALOGY 

In this section, we use the methods discussed in Ref. 
M to derive the hydrodynamic free energy of arbitrarly 
curved vortex lines. As usual, the vortex lines are viewed 
as imaginary-time world lines of 2D bosons. The need for 
a relativistic description is immediately apparent by con- 
sidering the self-energy part of the vortex energy, given 



in Eq. (1.2), where the arclength of a curve segment can 
be parametrized as ds = *J (dz) 2 + (dr±) 2 . Eq. (1.2) 
is formally identical to the imaginary-time action of a 
relativistic free particle of rest mass m <-» ei |32], 



S 



ds 



mc / V / c 2 (dt) 2 - {dr ± ) 2 , (2.1) 



where s\ and S2 are two space-time events and cr <-> z, 
with t = it the imaginary time. The speed of light c does 
not have a counterpart in the flux-line free energy, but is 
explicitly kept in this section to enable us to obtain the 
nonrelativistic limit. The correspondence between boson 
and vortex variables is summarized in Table 1. One im- 
portant diference between our mapping and the mapping 
to 2D nonrelativistic bosons discussed extensively in the 
literature is that here the boson mass, m&, corresponds 
to the vortex line energy rather than to the titl energy 
per unit length. This distinction becomes important in 
the case of anisotropic materials, as discussed in the Ap- 
pendix. 



TABLE I. Correspondence between boson and vortex vari- 
ables. 



bosons 


vortices 


CT 


z 


f3hc 


L 


s/h 


T/k B T 


rn.bC 2 /he 


ei/fcsT 


el /he 


4neo/kBT 


n/hc 


H 4>o/(4irkBT) 


Xb 


A 


U/c 


t « B/0o 
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As shown below, the vortex lines can be interpreted as 
the world lines of relativistic spin-0, charged bosons in 
2 + 1 space-time dimensions, interacting via a screened 
electromagnetic interaction p3]j . By applying standard 
methods |34||55| of rewriting the boson Hamiltonian in 
the language of second quantization and transforming to 
coherent states, the imaginary-time path integral repre- 
sentation for the grand-canonical partition function of 



5 r [$,$*,a,A] = J dr J d 2 r ± 



rribcf 



2 m h 



2 



1 



the bosons is written as 



^bos 



S»(r_ L ,/3B)=**(i-_ L ,0) 
#(r_ L ,/3B)=*(r_ L ,0) 



with the action 



ax)*| : 



2mfcC 2 
- CpW, A 



[(—ihd T + ebQo + ip r )&] [{iTid T + eba + ip r )§* 



(2.2) 



(2.3) 



where 



C F [a,A] = ^{ (V x a) 2 + ^(V x a) • A + (V x A 



Stt 



(2.4) 



is the free field Lagrangian density and V = (~d T , Vj_). 
Here, $(r±,t) is a complex scalar matter field that de- 
scribes bosons with positive and negative charge, with e& 
the boson charge. The world lines of bosons of opposite 
charge, which are antiparticles to each other, correspond 
to flux-lines with opposite vorticity. The constraint on 
the fields $ and <j>* comes from the permutation sym- 
metry requirement and reflects boson statistics. The bo- 
son field $ is coupled to a massive electromagnetic field 
a = (do, ax) an d the gauge field A = (Aq, Aj_) pro vide s 
the s cree ning. Integrating out the field A in Eqs. (2.2) 
and (2.4) under the gauge V • a = gives a mass term 
|a| 2 / (87rA^). Under Lorentz transformations, the field <& 
is a scalar and the boson action has full Lorentz invari- 
ance. The action also respects a global U(l) symmetry, 
<f> — > $' = $e M , which leads to the conserved current, 
j = (iojx), given by 

j = [$*( - ihV + — a)# + $(thV + —a)**] , (2.5) 



2rrib 



S r [p,0,a, A] 



lit, 



dr 



-C f [sl, A] 



rribC 



8rribp 



(Vp) 2 



with V • j =0. The temporal component jo/c of the 
2+1-current is the charge density, which is the appro- 
priate conserved quantity in a relativistic theory. Fi- 
nally, we have introduced a chemical potential, p, r , which 
couples to the conserved charge density. This should 
be contrasted to the nonrelativistic boson model which 
maps onto a liquid of directed vortex lines and where 
the relevant conserved quantity is particle number. The 
quantum relativistic model naturally incorporates spon- 
taneous creation and annihilation of particlc-antiparticle 
pairs. The flux-line analogue of this is the creation of ori- 
ented vortex loops. It will become apparent below that 
the chemical potential corresponds to the external field 
H. 

In order to highlight the connection with the nonrel- 
ativistic boson model used by Nelson and collaborators, 
we rewrite the boson field in terms of an amplitude and 
a phase, $(r±,r) = p(r±_ , t) exp [i0(r± , t)] , with the 
result 



2mfc 



(ftVifl + —ax) 



IwibC 



■(hd T 8 + Cba + i(j, r ) 



(2.6) 



A boson condensate phase is signalled by a macroscopic 
occupation p c of the p = momentum state, defined in 
terms of the average order parameter as 



p c = |(*(rx,r))| 2 



(2.7) 



Notice that in general p c ^ p = (p(r±,T)) = 
(|$(r_|_, t)| 2 ), although the two quantities are equal 
within the mean field approximation described below. 

In the spirit of Landau-Ginzburg mean field theory, we 
evaluate the partition function by the method of steepest 



descent. The stationarity condition gives nonlinear dif- 
ferential equations for the various fields. Restricting our- 
selves to solutions which are stationary and spatially ho- 
mogeneous, the extrema conditions yield the equations, 



SS r 
1p~ 
5S r 
Sa± 



TOftC 



rribC z 



2to{,c 2 
ax 



4ttA 2 



2mc 2 
= 0, 



(e b a Q + ip r ) 2 



0, (2.8) 
(2.9) 
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SS r e b 

t — = ?p(eba + i(i r 

oao rribc z 



4ttA 2 



0. 



1 / OF c \ 

(2.10) ^/c)^--(^) Tn = ^e bP0 , (2.13) 



One finds of course aj^ = 0, corresponding to the absence 
of charge current in the homogeneous state. A spatially 
homogeneous saddle point solution with |("I>}| 2 = po ^ 
exists only for \p r \ > nibC 2 j36|. The solution is given by 

\p r \-m b c 2 2 
tor \(j, r \ > rn.bC , 



Po 



4v\ 2 e 2 



Po = 0, for \p r \ < m b c 2 . (2-H) 

Within this mean field approximation, the occupation of 
the zero momentum state is p c = pa. The free energy in 
the condensate phase of bosons is given by 

(m b c 2 - p r ) 2 



F£(T,n,p r ) = ~k B T\YiZ c G = -hfi- 



8n\ 2 e 2 



(2.12) 



where f2 is the area of the 2D boson gas. The mean 
charge density, (jo/c), in this mean- field approximation 
is 



5 r [p,6»,a,A] = J dr J d 2 r± 



\ %vnbCS 



where the positive (negative) sign should be chosen for 
p r > (p r < 0) . For charged bosons the chemical poten- 
tial is associated with the electric charge and allows for 
a positive (or negative) average charge density. Antipar- 
ticles have a charge and a chemical potential opposite to 
that of particles. 

Within mean field theory, all the charge is in the con- 
densate phase. When fluctuations are incorporated in 
the theory, the mean charge density, (p) , differs from the 
condensate fraction, p c , even at zero temperature. 

The nonrelativistic limit for the action can be obtained 
by shifting the the chemical potential by to^c 2 , 



p r = p + VTlbC , 

where p is the nonrelativistic chemical potential |3£ 
action becomes 



(2.14) 



The 



-{drpf + z ( V -LP) 2 + ip{hd T 6 + e b a + ip) 

p 8m b p 



.^_(fiV_L0 + -a±) 2 + 7- ^(hd T 6 + e b a + ip) 2 + C f [sl, A] 
Zrrib c Zm b c z 



(2.15) 



The nonrelativistic boson model of Tauber and Nelson the 2D bosons only interact with the scalar field ao, and 
p5| is then recoverd by letting c — > oo. In this limit there is no magnetic interaction, 



f l3h f ( ft 2 P 
S nr [p,9] = dr d 2 r±< ihpd T 8 + ipe b a - pp + (Vj_p) 2 + - — (Wj 

Jo J { 8m bP 2r "fc 



il^Vim^-^l. (2.16) 



8tt lv ' J 87rA fc 



The scalar potential ao mediates the instantaneous 
screened Coulomb interaction. By integrating it out we 
recover the familiar instantaneous interaction term with 
Fourier components |/9(q)| 2 47rA 2 e 2 /2(l + q^A 2 ). 

In order to obtain the hydrodynamic free energy, 



we proceed as in Ref. [|| and perform a Hubbard- 
Stratonovich transformation to eliminate the boson am- 
plitude and phase fields in favor of familiar h ydro dynamic 
fields. First we eliminate the phase 6 in Eq. (2^) in favor 
of a new vector field P, with the result 



<S;[p,P,a,A] = / dr d 2 r 



lit, 



+i 



PPz 

rribC 2 



rribC 



I \ le bPr, r A 1 

(ebdo + ip r ) H P± ■ sl± + £f[&, AJ 



8rribC 2 P 



(drP? 



mbC 



8m,b P 



(V^p) 2 



m b po 



2m b 



p2 
_ z 
n 2 



(2.17) 



with the constraint density of the theory, which in turn is related to the tern- 

1 



c 



poral variations of the phase. The last term in Eq. (2.17) 
2 d T (pP z ) + Vj_ • (pPjJ = . (2-18) arises from the Jacobian of the functional integration and 

represents the nonlinear "ideal gas" part of the free en- 
Notice that the field P has three components. The tern- ergy. Finally, to make contact with the hydrodynamic 
poral component P z is related to the conserved charge fields used in our earlier work, we let 
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^ _ 29) an d integrate out the fields a and A to obtain an effective 

m b c ' ' action 



w.>± 



PF± (2.20) 



m b 



7 [2 %m b c A p srribp 2p c 

1 \ 47rA?e?/c 2 . . . , 

-2p^TW |t(q)l (2 ' 21) 



with the constraint 



V-t b = 0, 



(2.22) 



where q = {q z /c, q_i_) . 

Since we are interested in the fluctuating vortices 
rather than in the relativistic bosons, we shall rewrite 



ei 



{k B Tf 
8e lP 



(Vp) 2 



2p 



Eq. (2.21) in the vortex language using the "transla- 
tion" Table 1. We are interested in the thermodynamic 
limit where periodic or free boundary conditions give the 



same result and J Q cdr J d 2 r± — » 
L is the thickness of the sample, 
free energy is given by 



H, 



090 A 



4tt 



tz > + 



1 



47re A 2 



2VLL ^ 1 
q 



|t(q) 



St = So S dr ±> Where 

The resulting vortex 



(2.23) 



with the constraint given by Eq. ( 2.22 ) above. The fluc- 
tuating field t corresponds to the coarse-grained vortex 
density which determines the magnetic induction B in 
the material. The chemical potential, fi r , of the rela- 
tivistic bosons corresponds to the external applied field 
which gives rise to a net density of field-induced lines. 
The physical interpretation becomes clear by discussing 
the mean-field saddle point solution, obtained by apply- 
ing the stationarity condition to the free energy, with the 
result 











t=to,/3=PO 



= =>■ t z0 = ±cp , 



t=t ,p=p 



= 



H = ±(H, 



0o Po 



(2.24) 
(2.25) 
(2.26) 



\ / t=t ,p=/0o 

Conversely, the equilibrium solution is given by 

\H Q \ - H c i 

Po = , for \H Q \ > H c i, 

4>o 

Po = 0, for \H \<H cl , (2-27) 

where H c iei4ir / '<fto <-> mc 2 . In other words, po is a mea- 
sure of the number density of directed field-induced vor- 
tices. The field p is a scalar and it is defined to be always 
positive (for simplicity, in the following we will consider 
the case of H > 0, corresponding to the + signs in 



Eqs. (2.25) and (2.26). The field t z is proportional to 
the z component of the magnetic induction. Its equilib- 
rium value in a spatially homogeneous system is simply 



proportional to po because on large scales the contribu- 
tion from vortices precisely cancels that of antivortices, 
and t z o = cpo = cB z q/4>o . Locally, the fields t z and p 
can, however, fluctuate independently, allowing for spon- 
taneous vortex loop fluctuations, independent of the ex- 
ternally applied field. As will be seen more explicitely 
in the next section, the field p mediates the renormaliza- 
tion of the single- vortex stiffness due to such spontaneous 
loop fluctuations. 

The hydrodynamic free energy given in Eq. ( 2.25 ) 



provides a starting point for describing the long wave- 
length properties of a liquid of interacting {directed) field- 
induced vortex lines and oriented vortex loops, ft should 
be stressed that the distinction between directed lines 
and loops, while physically appealing, is strictly a single- 
line notion and loses much of its meaning in a contin- 
uum theory. At the level of hydrodynamics, spontaneous 
vortex loop fluctuations are incorporated via the field 
t% = t z — p. As we will see below, one can construct 
an effective theory where loops are integrated out. Their 
role then enters as a renormalization of the vortex-line 
tension. 



III. CORRELATIONS 

It is useful to evaluate the two-point correlation func- 
tion of the hydrodynamic fields appearing in the free en- 
ergy of Eq. (2.22). This is easily done within a Gaussian 
approximation for the free energy, obtained by introduc- 
ing fluctuations of the fields about their equilibrium val- 
ues, <5t = t — to and 5p — p — po, and expanding the 



6 



free energy to quadratic order in these fluctuations. For po > 0, the Gaussian free energy is given by 



??[6t,6p}=F 



2QL ^ 1 Un 



-Po 



Po 



~A 1 



l^(q)l 5 



^ + v( q )]\st( q )\ 2 

Po J 



[<^(q)fy(-q)+<^(-q)<5p(q)] , 



(3.1) 



where averages over the free energy are to be evalu- 
ated with the constraint 



Tf[St z ,t ± ] 



£q-£t(q) = 
1 



2flL 



E 



(3.2) 

(fc B T) 2 g 2 /4p 
ei + {k B T) 2 q 2 /4ei 



In Eq. ( |3.l| ), fo is the equilibrium value, V(q) — 
47re A 2 /(l + g 2 A 2 ), and q 2 = q\ + q 2 . It is instructive 
to integrate out bp to obtain the effective free energy, 



V(q) 



l^(q)l 5 



V(g) 



tx(q) 



(3.3) 



which is written entirely in terms of fluctuations in the and transverse to q^L = qj_/<Z_L according to 
local induction, as at long wavelengths, 

t ± (q)=q x ii(q) + (ixq ± )tl(q), (3.5) 

ffl«*)(ti,^j ■ (3.4) 



and using the constraint (3.2) to eliminate t\ in favor of 
Finally, by separating tj_ in its components longitudinal 8t z we obtain, 



2ilL 



till 



ei (fc B T) 2 g 2 /4e 2 

poql %ol + (k B T) 2 q 2 /4e 2 



v( 4 



l^(q)l 



^ + n<z)liii(q)i 2 



(3.6) 



Comparing the effective free energy of Eq. (3.6) to the obtained neglecting relativistic effects and given by 
conventional Gaussian hydrodynamic free energy [po|j7| 



T G [8t z ,t ± ] = 



1 



2nL ^ 

q 



Poq 2 ± q± . 



I^(q)l 



e -L + V(q)]\tl( q )f 
Po J 



(3.7) 



it is evident that "relativistic" effects yield short wave- 
length corrections to the single-vortex effective tension. 
Similarly, in the boson system they are responsible for 
corrections to the quasi-particle spectrum due to spon- 



taneous particle-antiparticle pair creation, Again, these 
effects arc important only at finite wavevector. 

It is now straightforward to evaluate the Gaussian two- 
point correlation functions, with the result 



(*t,(q)Jt,(-q)>< 



p k B Tq 2 



PoV{q)q 2 + e-iq 2 z 



2 JfcsT)V/4e? 
e l«_L l + (fc B T)V/4^ 



(3.8) 



(Sp(q}Sp(-q)) G = 



p k B Tq 2 



1 + PoV{q)/e 1 



nnV(a)nl + flff 2 + Fl „ 2 (^V/^? 1 + {k B T) 2 q 2 /4e 2 ' 
PoV{q)q +ei9 2 +ei9±T+pr^T)V74^ 



(3.9) 



(St z (q)Sp(-q)) G = 



p k B Tq 2 



P0V(q)q 2 + e 1 q 2 + e 1 ql 1 ^f^ 1 + ( W, 2 /4 £ 2 ' 



(3.10) 



7 



The correlations of the in-plane part tj_ of the tilt field 
are given by 



(t J (-q)^(q)> G =r«(q)^(q ± )+T«(q)^(q ± ) 



with 



and 



^(q) 



Pok B T 
ei + poV(q) 



l2(q) = 4<*.(-q)t,(q)> G 



(3.11) 



(3.12) 



(3.13) 



-P^(q-L) = q±iQ±j and fy(qi) = - P^ (qj_) arc the 
in-plane longitudinal and transverse projection operators 
respectively. To Gaussian order, the transverse part of 
the tilt field autocorrelator is the same as that obtained 
from the familiar hydrodynamics of directed lines. 



The structure function, given by the autocorrelator of 
fluctuations in the vortex line density, St z , differs from its 
"non-relativistic" counterpart only at finite wave- vectors. 
In the hydrodynamic limit it simply reduces to the famil- 
iar result obtained using the ( "non-relativistic" ) hydro- 
dynamics of directed lines, 



(St z (q)St z (-q)} G 



p Q k B Tq 2 A 



PoV(q)q 2 



(3.14) 



When "relativistic effects" are included the two fields 
p and t z are no longer independent. It is then instruc- 
tive to also consider fluctuations in a new field defined 
as their difference as t z — t z — p. This field fluctuates 
about a zero equilibrium value and vanishes identically 
when relativistic effects are neglected. Fluctuations in t\ 
may be interpreted as a measure of fluctuations due to 
the spontaneous excitation of vortex loops. Its correla- 
tion function is given by 



(^(q)^(-q))c 



p k B T 



q 2 z +PoV(q)q 2 /ei 

PoV{q)q +eig z +eig ±1+(fcBT)2g2/4e i \ t> , h i i 



(3.15) 



The long wavelength limit of this correlation function 
is given by 



lim lim <$#(-q)<$#(q)) G 



k B T 



(3.16) 



irrespective of the order of the limits (q z — ► first or 
q±_ — > first). In other words we have identified a corre- 
lation function that at long wavelengths yields the inverse 
single-line tension, e\. Calculating perturbative correc- 
tions to this Gaussian correlator will be a way to probe 
the renormalization of the single-line tension. 



IV. DISCUSSION 

We have presented a long wavelength (hydrodynamic) 
description of a liquid of arbitrarly curved vortex line 
and loops that describes on the same footing both field- 
induced and spontaneously generated vortices. The long- 
wavelength model was obtained by exploiting the map- 
ping of such a vortex liquid onto a gas of relativis- 
tic charged bosons in 2D that was recently discussed 
by Tesanovic 10 ll[] and by Sudb0 and collaborators 



[18, [L6|fl7| . Although our model and that of Tesanovic 
[11 1 yield free energies that are formally very similar, the 
two models differ in one important respect. In Ref. Jll| , 
the matter field $ that describes the relativistic quan- 
tum bosons corresponds to a fictitious vortex system of 
zero total vorticity. The field-induced vortices are sep- 
arated out of the field In our model, in contrast, 



the boson field <& corresponds to the actual flux-line sys- 
tem, including both field-induced and spontaneously gen- 
erated vortices. The externally applied magnetic field H 
which yields a non-zero vorticity enters explicitely in the 
free energy and controls the net mean flux threading the 
superconductor. It corresponds to a chemical potential 
for the bosons which yields a non-zero average charge. 
Furthermore, Tesanovic's fictitious vortex loops have a 
short-range scalar steric repulsion which enters as a $ 4 
coupling in the action. This is in addition to the usual 
screened magnetic interaction and results from the sin- 
gular gauge transformation which gives rise to the ficti- 
tious vortex loops. Such a steric repulsion does not exist 
among the actual vortex loops Jl|,^2| and therefore it is 
absent in our model. In our action, vortices only interact 
through the long range magnetic interaction (screened 
by the Chern-Simons term) which is cutoff at distances 
< £ to ensure a finite repulsive energy barrier at short 
distances. 

An important property of the flux- line array which pro- 
vides a direct measure of vortex correlations along the 
direction of the applied external field is the tilt modulus. 
A drawback of the hydrodynamics of directed flux-lines 
familiar from the literature j^] is that it has been impos- 
sible to separate the renormalization of the single- vortex 
part of the tilt modulus (related to the single-line tension) 
from that of the compressional part. The hydrodynamic 
model that we present here allows us to probe the renor- 
malization of the single-line tension through the autocor- 
relator of the t% field. This is an interesting direction for 
future work as the vanishing of the effective line tension 
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is considered to be a signature of the ^-transition sug- 
gested by Tesanovic jll] and by Sudb0 and collaborators 

Bps 

Finally, the hydrodynamic model of directed flux-line 
liquids has successfully been used to evaluate the renor- 
malization of the tilt modulus from columnar defects par- 
allel to the external magnetic field p7| . This type of 
quenched disorder suppresses the wandering of the flux- 
lines away from their average direction and it renders 
the directed line approximation appropriate. In contrast, 
other types of quenched disorder, namely point defects or 
splayed columnar defects with large splay, enhance vor- 
tex wandering, so that the resulting vortex liquid can no 
longer be described as a liquid of directed lines. The 
hydrodynamic model introduced here is particularly ap- 
propriate for studying the effect of this kind of disorder, 
which can be modelled by a random potential coupled 
to the field p. This is another interesting direction for 
future work. 
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APPENDIX: ANISOTROPIC 
SUPERCONDUCTORS 

Our discussion has so far been limited to isotropic ma- 
terials. High-T c superconductors are, however, layered 
materials where the anisotropy can play an important 
role and change substantially vortex behavior. For this 
reason in this section we generalize our model to include 
finite anisotropy. As usual, anisotropy is incorporated in 
the Ginzburg-Landau free energy via an anisotropic ef- 
fective mass tensor which leads to different values for the 
penetration and coherence lengths in the ab plane and in 
the c direction 

We first consider a uniaxial superconductor in an ex- 
ternal field H applied along the c axis chosen as the z 



direction. The energy of a single vortex line with po- 
sition parametrized as r = [rj_(z), z] (r±(z) can be a 
multivalued function of z to allow for the possibility of 
overhangs) wandering between points a and b is given by 



U 



self ~ £l 



f 



(dz) 2 + -j(drx) 



(4.1) 



where p — X c /X a b is the familiar anisotropy parameter. 
In copper-oxide high-T c materials, p >> 1. By compar- 
ing Eq. (4.1) to the action of a free relativistic particle 
given in Eq. (2.1), we see that a vortex line in a uniaxial 
material can be interpreted as the world line of a relativis- 
tic particle which moves in an "anisotropic space-time" 
with rescaled spatial coordinates, rj_ — > r±/p . At the 
single-vortex level, anisotropy effectively enhances the 
"relativistic" behavior, as it reduces the energy per unit 
length associated with transverse vortex fluctuations, as 
expected for a layered material. In this context, it may 
be tempting to interpret anisotropy as responsible for an 
effective reduction of the speed of light in the relativis- 
tic 2D boson problem with c — ► c/p. This interpreta- 
tion is, however, misleading as it does not carry through 
when interactions are included. In the interaction part 
of the vortex-line free energy, anisotropy does not simply 
lead to a rescaling of the transverse coordinates. This 
is apparent from the fact that the collective part of the 
wavevector-dependent elastic constants of a vortex lat- 
tice in a uniaxial material is not simply obtained from 
the corresponding elastic constants of an isotropic mate- 
rial with the replacement q± — » qj_p . 

In the interaction and the free field part of the boson 
model, the role of anisotropy is simply that of allowing 
for different scalar and transverse interactions among the 
bosons, precisely as originally proposed by Feigel'man 
and collaborators p9j . Note that the aforementioned 
parts in our relativistic boson action are formally the 
same as the corresponding parts in Feigel'man's boson 
action - our model differs in the free particle part only. 
The relativistic boson action that maps onto the free en- 
ergy of interacting vortex lines and loops in a uniaxial 
superconductor with an external field applied along the 
c axis is given by 



[P n r f TOC 2 i 
5 r [$,$*,a, A] = J dr J rf 2 r ± | -^-<M>* + ^ g2 [{-iW T + e b a + i(J, r )$] [(ihd r + e b ao + 



~|(-ffiVi + ^a ± )$| 2 + Cp [a, A] 

ZlTlb C 



(4.2) 



where 



Cf fa, Al 



\z x (-d T a ± - VooYI 2 + t~(V x a) • A + (V x A)' 

C Ab 



(4.3) 
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We stress that in this mapping the boson mass, nib, corre- 
sponds to the vortex line energy, e%, precisely as indicated 
in Table 1. This is indeed the appropriate interpretation 
as it is made apparent by comparing the single-vortex 
energy given in Eq. (|4.l| ) to the action of a relativistic 
boson, given in Eq. ( |2.l| ). In contrast, in all the previ- 
ous literature, and particularly in the work by Nelson and 



coworkers, it is the tilt energy per unit length, t\ = ei/p 2 , 
that is interpreted as the boson mass. 

Finally, by using the methods described in Section III 
and the translation Table 1, one can immediately obtain 
the hydrodynamic free energy of a liquid of arbitrarly 
fluctuating vortex lines and loops in a uniaxial supercon- 
ductor with H || c. It is given by 



{k B Tf 
8e x p 



oeip 2p 2p 2 



2n ^ 

q 



47re Ai 



«2 \2 i „2 „2 \2 



|t±(q) 



(l + q*Xl)(l + qf\i + qip*\l 



l**(q) 



(4.4) 



with ii = e\/p 2 and the familiar constraint V • t = 
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